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We study the behavior of an infinite system of ordinary differ- 
ential equations modeling the dynamics of a metapopulation, a set 
of (discrete) populations subject to local catastrophes and connected 
via migration under a mean field rule; the local population dynamics 
follow a generalized logistic law. We find a threshold below which all 
the solutions tend to total extinction of the metapopulation, which is 
then the only equilibrium; above the threshold, there exists a unique 
equilibrium with positive population, which, under an additional as- 
i^h | sumption, is globally attractive. The proofs employ tools from the 

C$ . theories of Markov processes and of dynamical systems. 



1. Introduction. The simplest models of population growth and regu- 
lation are formulated in terms of a more or less isolated population in a 
single habitat. However, the importance of the spatial dimension has been 
recognized in a number of ecological processes, resulting in one of the most 
(S| ■ active topics in theoretical ecology: see, for instance, the two recent collec- 

]Q ! tions [26] and [8] and the review article by Neuhauser [22]. These ideas have 

^ | stimulated the development of spatially structured stochastic populations 

models, as in [23] and [10], whose mathematical analysis is generally very 
hard. 

' A very simple model recognizing the spatial dimension of ecological pro- 

cesses was introduced by Levins [19] in 1969. He envisaged a metapopulation 
consisting of many distinct habitat patches, within each of which the pop- 
ulation behaves much as in the single population models, but which are 
^ . linked to one another by migration. In his highly simplified model, patches 

H ' are designated as occupied or not, and all occupied patches are taken to 
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be equivalent, irrespective of the number of individuals present. With these 
simplifications, he obtained a single differential equation, 



describing the behavior of the system: here, p = p(t) represents the propor- 
tion of occupied patches, is the extinction rate and cl is the colonization 
rate per occupied patch. Hence, an equilibrium exists only if cl > vl, and, 
in that case, the proportion of empty patches at equilibrium is vl/cl- His 
ideas have been widely used, both in theoretical papers and in wildlife man- 
agement problems (see, e.g., [15]). 

Levins' metapopulation model has two major weaknesses: on the one 
hand, it is based on a mean field assumption (the colonization rate in a 
patch depends only on the overall proportion of patches occupied); on the 
other hand, all patches are assumed to be equal and described simply as 
empty or occupied, disregarding local population dynamics. Addressing the 
first issue requires the consideration of spatial stochastic processes as men- 
tioned above. For the second, some authors have generalized Levins' model 
by taking into account the numbers of individuals in the occupied patches, 
giving rise to the so-called structured metapopulation models [12]: they con- 
sist either of a finite [21] or infinite number of ordinary differential equations 
[5], or of a partial differential equation [12, 13], where the structuring vari- 
able x represents the number of individuals per patch. However, very few 
analytical results are available for models of complexity comparable to ours, 
and the behavior of these models has mainly been explored through simu- 
lation. 

In this paper we investigate the deterministic approximation to the metapop- 
ulation model discussed in [1]. This is a stochastic mean field metapopula- 
tion model, in which the number of individuals in a patch is governed by 
a birth, death and catastrophe process, with the same transition rates in 
each patch, together with migration between the patches with a uniform 
transition rate 7 per individual, destinations being chosen uniformly at ran- 
dom among all patches. This last, mean field assumption is probably the 
least biologically realistic, but has been used in several papers [15], and may 
make very good sense for metapopulations of parasites in which the patches 
represent host animals. At all events, it makes the mathematical treatment 
substantially simpler. 

As is shown in [1], when the number of patches becomes very large, one 
can approximate the stochastic model with the following infinite system of 
differential equations: 





^ = c L p(l-p) - u L p, 



00 



Pi(t) 



(pi + di + j)i + v + P7 XI 3Pj (*) Pi (*) 
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oo 



+ b i - 1 (i-l) + fr/Y^JPjit) Pi-l(t) 
3=0 



(1.2) 



+ [d i+1 +-f](i + l) Pi+1 (t) 



i>l 



oo 



oo 



Po(t) = v[ X>;(t)-j>o(t) +(di + 7)Pi(*)-P7 Po(t), 



in which Pi(t) denotes the proportion of patches that are occupied by i indi- 
viduals, i > 0. The parameters bi and di represent the per capita birth and 
death rates in a patch occupied by i individuals, the catastrophe rate is v in 
each patch, the migration rate is 7 per individual, and p is the probability of 
a migrant successfully reaching another patch. Note that this model is very 
similar to those studied by Metz and Gyllenberg [21] as structured metapop- 
ulation models with finite patch size, and by Casagrandi and Gatto [5]. 
We also assume the following: 

(HI) ibi is concave and nondecreasing; idi is convex and nondecreasing. 

It can easily be seen that (HI) implies that bi is nonincreasing and di non- 
decreasing. Hence, there exist b^ = lim^oo bi and doo = lim^oo di, for which 
we further assume that 

(H2) b oc <d O0 + j(l- p) + v. 

Generally, in logistic demography, the existence of a carrying capacity is 
assumed: that is, there is a value K such that bx = dx, which automatically 
implies that 600 < doo- (H2) is weaker than that, and is, in fact, the natural 
condition: if > d^ + 7(1 — p) + u, there can be no nontrivial equilibrium, 
as is proved in Proposition 3.4. 

The assumptions of concavity of ibi and convexity of idi are satisfied in 
many examples, but not in all; for instance, a Ricker-type birth function 
bi = &oexp{— Pi} is not allowed. However, they are mathematically conve- 
nient assumptions, if the uniqueness of any nontrivial equilibrium solution 
to equations (1.2) is to be guaranteed, and we make use of them in several 
steps of our proofs; they could certainly be relaxed, but it is not easy to see 
what general conditions would better replace them. 

The existence and uniqueness of the solutions to (1.2) are established in 
[1], and a summary of those of her results relevant to this paper is given 
at the end of the section. In this paper we consider the possible equilib- 
ria 7r of (1.2), using stochastic coupling arguments that are developed in 
Section 2. There is always the "extinction" equilibrium, with 7r(0) = 1 and 
ir(i) =0, i > 1; this is also the eventual limit of all finite patch stochastic 



3=0 



3=0 



p(0) 
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systems, and makes the theory of quasi-equilibria of essential importance 
for such models. In addition, if a threshold condition is satisfied, we show 
in Section 3 that there is a unique nonnegative equilibrium having 7r(0) < 1 
(Theorem 3.1). This distribution is shown to be the equilibrium distribution 
for the single patch dynamics, in which immigration from outside is fixed 
at a constant "effective" rate, determined by the nonzero solution of a fixed 
point equation (3.3). In Theorem 4.5 of Section 4, we prove global conver- 
gence to this equilibrium when the threshold condition is satisfied, under the 
additional assumption that doo < +oo. The proof of convergence requires a 
lemma (Lemma 4.2) which is of some difficulty, because the system (1.2) is 
infinite dimensional. Its proof is the content of Section 5. 

The results of our paper give a rather complete description of this infi- 
nite system (1.2) of ordinary differential equations. Similar problems have 
recently been studied in other contexts, such as coagulation-fragmentation 
equations [2, 18], although for systems of equations of rather different struc- 
ture. It is possible that our methods could be useful in other contexts as 
well. 

We conclude the introduction by outlining the results that we need from 
[1]. First, note that the system (1.2) can be written in a more compact way 
as 



(1.3) 



where 



Pi = + m)p% + Ai_iPi_i + m + ipi + i 

+ P7I J^jPi (P*- 1 -P*) + u [dioYlPj -Pi)- 

\j = ) V 7=0 ' 



p_i(t):=0 for all t, /io:=0 and Ao:=0; 
A, := bit and /jj := (dj + j)i for all i > 1. 

It is proved in [1] that (1.3) is a well-posed problem in the space m 1 defined 
by 



to 1 = \x = (x ,xi,...) 



T el 1 , Y / j\x j \<oo\, 



equipped with the norm 

\\x\\ m = \x \ + 



oo 

i\xi | . 
i=0 
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More precisely, if Q is the infinite matrix 

' hi, if i + 1 = j > 0, 

-{(bi + (di+j))i + u), tfi = j, 

(Q)ij = Qij = \ (di + if i ~ 1 = j > 0, 

f(l-*io) + (di+7)*ti, if j = 0, 

k 0, otherwise, 

and g.,- = — qjj, we define the operator A by 

2?(v4) = |u G to 1 :^g fc |n fc | < ooj; 

(1.4) 

Au = uQ; (Au)i = ^2q ki Uk. 

k 

Then it turns out that the closure A of A is the generator of a C -semigroup 
on m 1 (see also [24]). 

We then define the map F : m 1 — > m l by 

(1-5) F(p)= P1 ( f>,- Wxfp) -/(£)), 

\i=o / 

where (T_i(p))j :=pj_i and J is the identity. .F is Lipschitz and, in this 
notation, the system (1.3) can be written as 

/ p' = Ap + F(p); 

1-6 ~ ~ 

p(0)=p°. 

The following theorem is proved in [1]. 

Theorem A. For every p° > G f(-A) and any T > 0, iaere exists a 
unique p[t) > G C([0, T];T>(A)) n C 1 ([0, T]; to 1 ) satisfying (1.6). C/ear/y, 
p(i) u>iZ/ a/so satisfy (1.3) componentwise. 

Moreover, if p° G C = {p G to 1 :p > 0, Yl^aPj = 1}; ^en p(t) G C /or a// 
t>0. 

Since the Pi(t) represent the frequencies of sites with i individuals, the 
condition p(i) G C is quite natural, and most of the following results relate 
only to that case. 



2. Immigration, birth, death and catastrophe processes. The analysis of 
the differential equations system (1.2) is accomplished indirectly, using the 
properties of a number of associated birth and death processes. We make 
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several comparisons based on couplings of such processes, which exploit the 
fact that birth and death processes cannot cross without meeting. A good 
general reference is [20]; in particular, see pages 3 and 4. We begin with a 
simple lemma. 

Lemma 2.1. Fix a positive integer J, and let V = (Vt, t>0) be the birth 
and death process on the integers j > J with transition rates 

+ l at rate j<j>, j > J; 

(2.1) 

J ~>j ~ 1 at rate jfi, j>J+l, 

for some (j),fj, > 0. Then, ifE^ denotes expectation conditional on Vj) = m, 
for any j' > J , we have the following: 

1. Ifcj> < then E^{V 2 } < {jV/(At " ■ 

2. If(f)>n, t/ien E^"){y t 2 } < {2j>e^-^V(<^-^)} 2 . 

3. If <f> = n, then, for any e > 0, E^{V t 2 } < {2/(0 + e)e £t /e} 2 . 

Proof. It is enough to conduct the proof for / = J: for / > J, the 
F-process is stochastically smaller than a 1^-process defined with J replaced 
by f. 

Suppose first that (j> < fx, in which case V is positive recurrent. Observe 
that a monotone coupling of two realizations of ^-processes, one with initial 
state J and the other starting with its equilibrium distribution tt, shows that 
E( j )(V^ 2 ) < W{Vq) for all t. Now tt satisfies the detailed balance equation 

j</m(j) = (j + l)fj,ic(j + 1), j>J; 
hence, jvr(j) < J(c/)/[iy~ J for all j > J, from which it follows that 

e( j \v 2 ) <e*(v 2 ) <jj2jW»y- J 

= JMv - <f>)~ 2 + J*i*(j* - 

this proves part 1. 
If <j) > fi, we have 

(2.2) EV t = E(V t I[T 1 < t])+E(V t I[n > t\), 

where 

Ti = inf < t > 0; V t = J, max V s > J + 1 \ < oo. 
L o<s<t j 

Note that P( J )[n < oo] = fj,/<f> and that E^iVflln > t}) < E( J )(V; 2 ), where 
V is a birth and death process on Z + with rates as in (2.1), but now for all 
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j > 0; this latter bound implies that 

E^{V t 2 I[n > t}) < {E^V t } 2 + Var( J ) V t 

(2.3) <{Je^ )t } 2 + JVar^Vt 

< J 2 e 2 ^-^ + j e 2 (^) i {2(0 + n)/{<l> - fj,) - 1}. 

Also, again by a monotone coupling of two ^-processes, 

E^iVflin < t]) < P (J) [n < t]E( J \v t 2 ) = (fi/(f))E^ J \V 2 ), 
and hence, from (2.2) and (2.3), 

E (J) (V 2 ) < {0/(0 - A*)} J 2 e 2 ^{2{<P + AM - 
proving part 2, and also, once more by stochastic comparison, part 3. □ 

Now let X := (Xt, t > 0) be an immigration, birth and death process 
with per capita birth and death rates (3j and 5j, respectively, j > 1, and 
with immigration rate A. Suppose that the function n(3 n is concave and 
increasing in n > 0, and that n5 n is convex and increasing. Then it follows, 
in particular, that (3 n is decreasing and 5 n is increasing in n > 1; we define 

(2.4) c:= lim (3 n - lim <5 n . 

n-^oo n— >oo 

Theorem 2.2. Lei X and c be as above. Then: 

1. There exist constants C\ and Ci(e), e > 0, such that 

dil+j 2 ), ifc<0, 
Ci(e)e 2 ( c+£ )*(l+j 2 ), foranys>0, ifc>0. 

2. There exist constants C2 and 02(e), e > 0, such that, for all m > 0, 

C 2 , ifc<0, 
C 2 (e)e( c+£ )' , for any e>0, if c > 0. 



E^(X 2 ) < 



Q<lfc(™+l) Xt _ E (m) Xt < 



3. In either case, for all m > 0, 

E (m+l) Xt _ E (m) Xt > E (m+2) Xt _ E (m+1) ^ 

Proof. Let /3'- := /3j + j _1 A for j > 1. Then note that, for any positive 
integer J, a simple monotone coupling of two birth and death processes 
shows that, if Xq < J, then X is stochastically smaller than a birth and 
death process V as in Lemma 2.1, having (j) = (3'j and /x = 5j and starting 
with Vq = J, since Vq > Xq and the sequences /3j and <5j are nonincreasing 



<s 
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and nondecreasing in j, respectively. If c < 0, choose J so that f3'j < Sj, and 
use Lemma 2.1 part 1 to give 

E V)(X?) < ^ J \v 2 ) < {J8j/(5j - j < J, 

(2.5) 

E^(X 2 ) < E®(V?) < {jSj/(Sj - P'j)} 2 , j > J. 
If c > 0, choose J so that 

5j < 0'j < 5j + c + e, 

if this can be done, and use Lemma 2.1 part 2 as above to give 

(2.6) E^(X 2 ) < E^(V 2 ) < {2m a x{J,j}(3' J e^ t /((3' J - 5j)} 2 . 

The only remaining case occurs when A = and the sequences (3j and 5j are 
both constant for all j > J for some J, in which case Lemma 2.1 part 2 or 3 
can be applied directly. Combining this observation with (2.5) and (2.6), 
part 1 is proved. 

We now turn to part 2. Let (Y, W) := ((1^, Wt), t > 0) be a two-dimensional 
pure jump Markov process with transition rates given by 







at rate 


ifr + A, 








at rate 


iSi, 




(iJ)- 




at rate 










at rate 


(i + j)8 i+j 





for all i,j > 0. All the transition rates are nonnegative, because both nf3 n 
and nS n are increasing. Then the processes Y and Y + W are also Marko- 
vian, both having the same generator as the immigration, birth and death 
process X. Thus, we can couple realizations X 1 and X 2 of X with Xq = m 
and Xq = m + 1 by realizing (Y, W) with Yq = m and Wo = 1, and setting 
X 1 = Y and X 2 = Y + W. Then it is immediate that X? - X} = W t > for 
all t; the next step is to bound EW*. 

However, just as before, a simple monotone coupling shows that W is 
stochastically smaller than a birth and death process V as in Lemma 2.1, 
having (j> = and fj, = Sj and starting with Vq = J, since Vo > Wo and, for 
any i > 0, 

(i + m +j - ifr < jfy < jft, j > j, 

and 

(i + j)<5j+j - i5i > jSj > jSj, j>J + l, 

by the concavity of n(3 n and the convexity of n5 n . Thus, in particular, EVFj < 
E^Vt, and the bounds on ~E( J \V 2 ) obtained in part 1 can be invoked, 
completing the proof of part 2. 
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For part 3, we extend (Y, W) to a four-dimensional pure jump Markov 
process ((Y, Wt, Ut, Vt), t>0) with transition rates 



n - 


+ n + e « 


at rate 


ifii + A, 




n - 


+ n-6« 


at rate 






n - 


+ nW 2) 


at rate 


(i + j)Pi+j 




n - 


^n- e ( 2 ) 


at rate 


(i + j)5 i+j - i8 if 




n - 


+ nW 3) 


at rate 


(i + k)P i+k -iPi, 




n - 


- n - e ( 3 ) 


at rate 


(i + k)5 i+k - iSi, 




n - 


^nW 4) 


at rate 


(i + 3 + - 


(i + j)P i+j 


n - 




at rate 


(i + j + l)S i+j+ i - 


(i+j)S i+j , 



when n = k, I) is such that k^l, the last four transitions being replaced 
by 

n -> n + e (3) + e (4) at rate (i + j + k) Pi+j+k - (i + j) Pi+j , 
n — > n — e^ 3 ) — e^ 4 ) at rate (i + k)5i + k — iSi, 

n^n + e (3) at rate (i + k)P i+k - ifo - (i + j + k)P i+j + (i + j)Pi+j, 
n^n-e (4) at rate (i + j + k)S i+j+k - (i + j)S i+j - (i + k)5 i+k + i6 i} 

when n = (i,j,k,k), all transition rates being nonnegative because of the 
assumptions on nb n and nd n ; here, e^ m ' ) denotes the mth coordinate vector. 
The four processes Y , Y + W, Y + U and Y + W + V are Markov, and each 
has the same generator as the immigration, birth and death process X. Thus, 
realizations X 1 , X 2 , X 3 and X A of X with Xq =m, Xq = Xq = m + 1 and 
X$ = m + 2 can be obtained from (Y, W, U, V) by setting X 1 = Y, X 2 = Y + 
U, X 3 = Y + W and X 4 = Y + W + V and taking Y = m ,W = Uo = V = 1. 
Thus, E< m+1 )Xi - E m A^ = EC/* and E( m+2 )X t - E^ m+ ^X t = EV t . Initially, 
Uq = Vq = 1. Thereafter, both ?7 and V make only unit jumps, and at any 
time at which U and V are equal, either they can jump together, or U can 
increase by 1 or V can decrease by 1. Thus, U is always greater than or 
equal to V, and, for each t > 0, Ut > Vt with positive probability. Hence, for 
all t > 0, 

E (m+l) Xt _ E m Xt = EUf > EVf = E (m+2) Xt _ E (m+1) x ^ 
proving part 3. □ 

The theorem above is used in the study of our main object of interest, 
a family of immigration, birth, death and catastrophe processes , indexed 
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by an immigration parameter s. The pure jump Markov process has 
transition rates 

j -»• j + 1 at rate q jtj+1 := jbj + p-/s, 

(2.8) j ->• j - 1 at rate e^-i := + 7)) 

j — > at rate q^o := 

and n6 n is assumed to be increasing and concave, nd n to be increasing and 
convex. The process starting with any initial distribution ip can be con- 
structed as follows from a sequence of independent realizations X^ , , . . . 
of an X-process with parameters f3i = bi, A = pjs and 5i = di + 7, and with 
X^ ~ if) and Xq 1 ^ = 0, n > 1. Let the times (T n , n > 1) of the catastrophes 
be the partial sums of independent negative exponentially distributed ran- 
dom variables (E n , n > 1) with mean 1/V, which are also independent of 
(iW, n > 0). Set N(t) : = min{n > 0:T n < i}, where T := 0; then define 

(2.9) z[ s) :=X^\t-T N{t) ). 

A pair of -processes Z^' 1 ^ and Z( s ' 2 ) with different initial states k > / 
can then always be coupled by using the same sequence of X-processes 
(X^ , n > 1) and taking XW = Y + W, X(°> 2 ) = Y, where (Y, W) is as in 
(2.7) and has Yq = I and Wo = k — I. With this construction, it is clear that 
z (s,i) > z {s,2) for aU ^ that p[ Z («,i) > z («,2)j < g-^t and that 

(2.10) < E(Z t (s>1) - Z^' 2) ) = e- ut EW t = e~ ut (E^X t - E®X t ). 
Defining 

(2.11) c := lim b n — lim d n — 7, 

n— >oo n— >oo 

and assuming that c < z^, it thus follows from Theorem 2.2 part 2 that, for 
k> I, 

< E (fc) z (s)_ E (Z) z (s) 

(2.12) Uk-l)C 2 e- vt , ifc<0, 
~ 1 {k-l)C 2 (l(v-c))exp{-l(v-c)t}, if c>0. 

Thus, if — > K is Lipschitz with constant K(f), then 

(2.13) \E^f(Z ( t s) ) - E^f(Z { t s) )\ < C{k - l)K{f)e~ at 

for some C, a > 0. Furthermore, from (2.10) and Theorem 2.2 part 3, we 
have 

(2.14) E^ m+1 ^z[ s) - E^Z\ S) > E( m + 2 >Z t {s) - E^ m+1 ^z\ s \ 
for all m, t > 0. 
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Theorem 2.3. Let be as defined in (2.8), with nb n increasing and 
concave, nd n increasing and convex. Suppose that c<u, where c is as defined 
in (2.11). Then, for s > 0, is positive recurrent, and its equilibrium 

distribution 7r^ has finite mean equal to lim^ooE^ ) z[ s ^ ; furthermore, for 
any < 5 < 1 for which c(l + 5) < v , we can find K\{5) < oo such that 

(2.15) E^{(4 s) ) (1+5) } < Ki(S){l + j {1+5] }, 

for all t>0 and j > 0. 

If s = 0, the state is absorbing for Z^ s \ and the only stationary distri- 
bution is 7r(°) = A{o}, giving probability one to 0. 

PROOF. The case s = is immediate, so we now suppose that s > 0. 
If v = and c < 0, the detailed balance equations 

(2.16) (jbj + frysfa = (j + l)(d j+1 + 7 )^ +1 , j > 0, 
are satisfied with 



bj + fyysj 1 
vtj+1 < — — ttj < (1 - e)irj, 



for some e > and for all j large enough, because c < 0. Hence, (2.16) have a 
nonnegative solution with geometrically decreasing tail, and the conclusion 
of the theorem follows. 

If v > 0, positive recurrence is immediate. Construct Z^ with Zq = 

from a sequence of X-processes as above. Then, if m(t) := m^ s \t) := Ez[ s \ 
we have the renewal equation 

m(t) = e~ ut E{X t \X = 0} + [* ue~ vu m{t - u) du. 

Jo 

Now, by Lemma 2.1 part 1, 

'Ci, ifc<0, 

C 2 ((z/-c)/2)exp{i(i/ + c)t} ; ifc>0, 

for suitable constants C\ and Oj. Furthermore, a monotone coupling of two 
X-processes with different initial conditions shows that E^Xt increases 
with t. Hence, the key renewal theorem ([11], page 363) can be applied to 
conclude that 

(2.17) m (s) (oo):= lim EZ ( t s) = v f°° e - vt E{X t \X = 0} dt 

t-*oo J 

exists and is finite. But now, because Z^ is nonnegative and positive re- 
current, it follows from (2.17) that 7P S ) has finite mean, satisfying 

^(e) =E w(s) (Z ( Q s) ) <m (s) (oo), 
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where e(j) := j for all j > and 7r«(/) := Efc^o^/C*)- 

Finally, for any < 5 < 1 for which c(l + <5) < i/, a similar renewal ar- 
gument can be employed, again appealing to Lemma 2.1 part 1, to show 
that 

m { s s \t) :=E(°){(Z t (s) ) (1+5) } 

(s) 

is uniformly bounded for all t; hence, the sequence of random variables Z\ 
is uniformly integrable, and thus, in fact, 

(2.18) 7rW(e) =mW(oo), 

proving the first two claims of the theorem. Noting also that, for any e > 0, 
K^{(Z { t s) f +5] } 

(2.19) = e~ l/t E^{X t (1+5) } + C ve- uu mf{t-u)du 

Jo 

< mj°(t) + (1 + j (1+5) )c[ 1+S)/2 (e) exp{(l + <5)(c + e)t - ut}, 
from Theorem 2.2 part 1, the remaining claim is also proved. □ 

With these preparations, we can now prove the main result of the section. 
The assumptions of Theorem 2.3 are still in force. 

Theorem 2.4. Let ir^ denote the equilibrium distribution of the pro- 
cess Z^; then ir^(f) is continuous in s for any Lipschitz function f . Fur- 
thermore, if e(j) := j for all j > 0, then 7P s '(e) is an increasing, strictly 
concave function of s . 

Proof. Let be the generator of the process Z^ s \ so that 
(A^h)(j) :=5>j{fc(0-Mj)} 

(2.20) = (jbj + Pis){h{j + 1) - h(j)} 

+ j(dj + l){h(j - 1) - h(j)} + u{h(0) " Kj)}, 

and, for any Lipschitz function / with constant K(f), let 6^ s \f) be defined 
by 

roc , . 

(2.21) e^(f)(j):=- {^f{Z { t s) )-^ s \f)}dt. 



We begin by showing that 9^ s \f) is a solution h to the equation 
(2.22) {A^h){j)=f{j)-^ s \f), j>0. 
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First, realizing ir^(f) = W is) f(z\ s) ), it follows from (2.13) that 
|EW/(Z t (s) )-7T (s) (/)l 

= E4 s) {E (i) /(4 s) )-^/(4 s) )} 

< CK(f)e~ at Y,4 S) \k-j\ < CK(f)e- at {m^(oc)+j}. 
Hence, 9^ 8 \f) given in (2.21) is well denned. Now set 

{ T\f)(j)--=- l_ {^ j) f(Z { t s) )-n^(f)}dt, 



o 



noting that lim T ^oo 0? ) (/) (J) = ^ s) (f)U) by (2.23). Conditioning on the 
first jump gives 



= "E (i) | [ T {f(Z { t s) )-7r^(f)}dt 



o 

-e-* T T{/(i)-vrW(/)} 

rT ( / f-T 



(2.24) 



I e " 9j "E^{ E( °(/ V(4 S) )-^ (S) (/)}^ 



= -q-\l-e^ T ){f(j)-^(f)} 
Now, by (2.23), 

for all T and it, and q~i > only for I = 0, j — 1, j : + 1. Thus, letting T — ► oo 
in (2.24) and using dominated convergence, it follows that 

/•oo 

{s \f)(j) = -qj l {f(j)-^ s) U)} + Y.^ / e _ * tt w (/)(*)<*«, 

#i 

or 

M - * {s) (f) = E^ (s) (/)(0 - ^ (s) (/)(i)i- 
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Thus, solves (2.22). 

Furthermore, again using (2.23), 

(2.25) e^(f)(j + 1) - e^(f)(j) = /°V%(4 S) ) - E^f(z ( t s) )}dt 

Jo 

and (2.13) immediately gives 

(2.26) |# (s) (/)(j + 1) - # (s) (/)(j)l < CK(f)/a; 

thus, A0( s \f) is bounded and, hence, also Lipschitz, with constant 

(2.27) K{A9^{f))<2CK{f)/a. 

Now, by Dynkin's formula ([14], Theorem 2), it follows that ^(A^h) = 
for all s, for any Lipschitz function h. In particular, for any t > —s, us- 
ing (2.20), 

(2.28) =E- (3+t) {(^)^ s )(/))(4 s) ) + /07 t(A^)(/))(Z«)} 

= E^ +i) {f(Z { s) ) - tt«(/) +p7i(A^(/))(4 S) )}- 
Thus, from (2.28) and (2.26), it follows that 

(2.29) \^ +t Xf)-^Xf)\< P1 t\\M^{f)\\< P1 \t\CKU)/a 

for any Lipschitz function /, so that 7P S )(/) is continuous in s, proving the 
first part of the theorem. 
It then also follows that 

| 7r (^)(/)_ 7r W( / ) + /97i7r W(A^ s )(/))| 

< pj\t\\^ s+t \A9^{f)) - TT^ s \A6^{f))\ 

and, hence, that 

(2.30) ^ W (/) = -P7T W (A0«(/)). 
Taking / = e, this last can be re-expressed using (2.21) as 

(2.31) ^ (s) (e) = PTlE 7 ^ / o °°{5(4 S) + l,t) 

where g(j, t) := K^Z^ . Hence, from Theorem 2.2 part 2, it follows that tt^ (e) 
is increasing in s, proving the next part of the theorem. 
Now, from (2.28) with u and 2u for t, 
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and 

7r(«+2«)(/) - ^(f) = -2 P1 u^ s+2u \Ae^(f)) 
giving, again from (2.28), 

= -2 P1 u{^ s+2u \A9^\f)) - ^ s+u \A9^ s \f))} 
= -2 P1 u{^ s) {A9 {s \f)) - 2 P1 u^ s+2u \A9 {s) {A9 {s) {f))) 
- 7T W {A9 ( - s) (/)) + P1 u^ s+U ^ {A9 {s) {A9 {s) (/)))} 

= 2( P1 u) 2 ^ s \Ae^ s \Ae (s \f))) + 7?, 

where 

M < lO(p7kl) 3 C^(A0^(A0W(/)))/a < 40(C P7 |n|/a) 3 K(/), 
this last by (2.29) and (2.27). Hence, ^ s \f) is twice differentiable, and 

(2.32) ^^(Z) = 2p 2 7V s) (A0«(A0« (/))). 

Now, using the formula given in (2.25), it follows that 
A0W(A0OO(/))( 



(2.33^ /*oo ^ /*oo 

= ^ |lE (m+1) I {g f (Z { t s) + l,w)-g f (Z { t s \w)}dw 

J^{g f {z[ s) + -^(Z t W ,«;)}d«;}dt, 

where g f (l,w) :=E^ f(zl s) ). To evaluate (2.33), realize Z^ with z{, s,1) = 

j + 1 and Z^ s ' 2 ^ with Z^ 2 ^ = j as before, using the Markov process (Y, W) 
of (2.7), with y = i and W = 1, so that 

Z( (»,i) = 2 t M + w[£l>t]i 

where E\ is an independent negative exponential random variable with 
mean Ijv. Thus, 

A9^(A9^(f))(j) 

(2.34) =E&) / / e- l '*{ ff/ (y t + Wt + l,t 1 ;)- 5/ (Y t + ^ t , U ;) 
jo jo 

- 5/ (y + 1, w) + 5/ (y, «;)} dw dt. 
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In order to use (2.32) to investigate the curvature of 7r( s )(e), we take / = e 
in (2.34). Then, for any k>l, 

g e (k + l,w) - g e (k, w) - g e {l + 1, w) + g e (l, w) 

= {^ k+l ^Z^ - E^4 S) } - {^ {l+l) Z^ - E®Z$] < 0, 

from (2.14), for all w > 0, so that the integrand is always negative. Hence, from 
(2.34), it follows that A(9W(A0( s )(e))(j) < for all j and s, and thus, 
from (2.32), 

^7r( s )(e)<0, s>0. 
This completes the proof of the theorem. □ 

3. Equilibria. We now investigate the equilibrium solutions of (1.3). For 
the sake of simplicity, we shall assume here and in all that follows that p = 1 . 
There is no real loss of generality in this, since one could set d\ = di + ^(l — p) 
and 7' = 7/3, and write (1.2) using d' and 7' in place of d and 7. In biological 
terms, unsuccessful migration is just one cause of death. 

If 7r € to 1 is such a solution, and 

00 
i=o 

then 7r must solve 

= - [(bi + d,i + ^)i + v + js]iTi 

(3.1) +[6;_i(i-l)+7s]7T i _i + [(i i+1 +7](i + l)7r i+ i, i>\- 
= 1/(1- 7T ) + (dl + 7)71*1 - 7S7T . 

Hence, tt must be the equilibrium distribution of the immigration, birth, 

(s) 

death and catastrophe process Z\ , which we studied in detail in Section 2. 
From Theorem 2.3, and using (HI) and (H2), we know that Z^ s > has a unique 
stationary distribution tt^ s \ which has finite mean denoted by 

(3.2) G ( S ):=n^(e)=J2j4 } - 

In order to have an equilibrium solution of (1.3), s must be equal to 7r^^(e); 
in other words, we look for a solution to the equation 

(3.3) s = G(s), 



a fixed point of the function G. 



A METAPOPULATION MODEL 



17 



Theorem 3.1. Suppose that (HI) and (H2) are satisfied. IfG'(0) > 1, 
then there exists a unique positive fixed point s* of G; if G'(0) < 1, then 
G(s) < s for all s>0. 

Remark 3.2. Note that s = is always a fixed point of G; the corre- 
sponding equilibrium distribution 7r(°> is the vector e° = (1, 0,0,... ) T , which 
can be interpreted as the extinction equilibrium. 

For the proof, we need a technical point. 

Lemma 3.3. Let 



I .7 



m 2 = |i£l 1 , j 2 \xj | < oo 



i/rai/i norm 



oo 

|x|| m 2 = |X | 



-2 1 

i=Q 



1 1 ■ 



and let be the part of A in m 2 , that is, 

T>(A 2 ) = {x£V(A):Ax€m 2 }; A 2 x = Ax. 
Then, if p(0) €T>(A2), p(t) satisfies 

oo 

(3.4) ]Tj 2 a>^) <oo. 

Proof. We first note that the restriction of e At to m 2 is again a C°- 
semigroup. This can be established following, with obvious changes, the 
proofs in [1]. In fact, repeating step by step the proof of Proposition 6.5 
of [1], one sees that A 2 — uj is dissipative, as long as u > 3maxj6j. The 
density of the range is then established exactly as in Proposition 6.6 of [1]. 

Moreover, repeating the proofs of Lemmas 6.8 and 6.9 of [1], one sees that 
the domain of the restricted semigroup is contained in the set 



jxe^ 1 , ^fd^x^ < ooj, 



and (3.4) follows. □ 



Proof of Theorem 3.1. First of all, it is proved in Theorem 2.4 that 
G is an increasing, strictly concave function in s > 0. We now establish two 
further properties of G: 

(3.5) G(0) = and lim G(s)/s < 1. 
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The first of these follows because, when s = 0, the equilibrium distribution 
is concentrated at 0, so that its mean is 0. 

For the limit as s — ► oo, we let m(t) = m^ s \t) = E(Z t ^), with z[ as de- 
fined in (2.8), noting that G(s) = lim^oo m( s \t) as shown in (2.18). Letting 

Pj(t) = ¥(Z^ = j), we can write m(t) = J2j JPj(t)- 

The forward equations satisfied by p(t) can be written as p'(t) = A s p(t), 
where 

A,:V(A)^m\ (A, p)i - ( <** + ^ ^ 

{ (Ap)o - 7spo, i = 0, 

is a bounded perturbation of the operator A defined in (1.4). 

Hence, if the initial value p(0) is in V(A), then p(t) =e Ast p(0) is differ- 
entiable function from R to m 1 and we have 

oo 

m'(t)=J2jPj(t) 
i=i 

oo 

(3.6) = £j{(0' " + 7«)Pi-i(*) 

i=i 

-((6i + (dj+7))i + 7« + ")Pi(*) 

+ (j + l)(d j - 1 +j)p j+1 (t)}. 

If p(0) € D(A2), the condition (3.4) allows the order of the sums in (3.6) to 
be interchanged, and, with some manipulations, we obtain 

(3.7) m'{t) =Y,j( b j ~ dj)Pj{t) " (7 + v)m{t) + js. 

j 

Using the concavity of xb(x) and the convexity of xd(x), we obtain 

(3.8) ^jbjPj(t)<m{t)b(m(t)) and ^jd j p j (t)>m{t)d{m(t)). 

j j 

Hence, from (3.7), it follows that 

(3.9) m'(t) < m(t) [b{m(t)) - d(m{t)) - 7 - v] + 7s, 

so that m(t) < x^ s \t), where x := x^ s '(t) is the solution of the Cauchy prob- 
lem 

x' = x[b(x) — d(x) — 7 — v\ + 7s, 

(3.10) 

x(0)=m(0). 

Since V(A2) is dense in T>(A), it follows that m(t) < x^ s \t) for all p(0) € 
V{A). 
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Set 

(3.11) a = v + — 6oo > 
because of (H2), and choose fh such that 

b(fh) — d(fh) — v = —a/2, 
if this is possible; otherwise, set fh = 0. In any case, we have 

(3.12) b(m) — d(m) — v < — a/2 for m > fh. 

Take s such that 7s = fh{% + 7). Then for all s > s, there exists such 
that x(t^) = fh and x(t) > fh for t > t^ s K Then, using (3.12), we have 

x'(t)<is-x{t)(\a + i) for t>r (s) . 

Hence, 

x(t) < me-W 2 ^^-^') + 7 s t e -(<V2+7)(*-) da 



7 S _ f ^ _ ^ e -(«/2+7)(t-rW) < ^ S 

7 + a/2 \7 + «/2 / -7 + 0/2' 



so that, using (2.18), 



7s 

G(s) = lim m(t) < lim x(t) < 



t^oo t^oo 7 + a/2 

and, hence, 

lim G(s)/a < ^— < 1, 

*^oo w/ 7 + a/2 

as stated above. 

Turning now to the fixed points of G, note that G(0) = and G is strictly 
concave; hence, G(s) = s has, at most, one other solution in s > 0. Since 
also lim^oo G(s)/s < 1, it follows that there is a unique positive solution of 
G(s)/s = 1 if G"(0) > 1; otherwise, if G'(0) < 1, we have G(s)/s < 1 for all 
s>0. □ 

The next result shows that assuming condition (H2) to be satisfied is not 
restrictive, when looking for positive equilibria of (1.2). 



Proposition 3.4. If (H2) is violated, there are no nontrivial equilib- 
rium solutions to (1.2). 
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Proof. If 7 = 0, the proposition follows immediately from the p -equation 
in (1.2). Otherwise, the process is stochastically larger than a pro- 
cess which has bj = and dj = doo for all j, and the same is true if 
bj = b 00 — b* for any < 6* < b^. Letting m := ¥}z[ s \ note that, as for (3.6), 

(3.13) m'(t) = m(t){&oo - b* - d^ - 7 - z/} + 7s = -m(t)(a + 7) + 7s, 

where a' = a + 6* and a is as in (3.11). Suppose now that a < 0, so that (H2) 
is violated. If a < —7, choose 6* so that a' = —57; otherwise, take 6* = 0. 
Then it follows from (3.13) that 

G(s) = lim EZ^ s) > lim m(t) = s{j/{j + a')} > s, 

t — >oo t — >oo 

for all s > 0, and there can be no s > for which G(s) = s. 

Finally, if a = and 7 > 0, then the condition c < v of Theorem 2.3 
is satisfied, so that, from Theorem 2.4, the function G is strictly concave 
and G(0) = 0. The argument above then gives G(s) > s for all s, which 
therefore precludes the existence of any s > with G(s) = s. This completes 
the proof. □ 

Remark 3.5. From (2.31), we see that 

(3.14) G"(0) = 7 r^Z^dt. 

Jo 

Thus, G"(0) is the average number of successful propagules produced in 
a patch colonized by a single immigrant, before population extinction in 
that patch, disregarding other colonizations. This number may be consid- 
ered a reproduction number for colonizers of an empty habitat, as used 
in epidemic models [9], thus, G"(0) > 1 is the natural condition to ensure 
(meta)population persistence. Indeed, a similar condition has been presented 
by Chesson [7] and Casagrandi and Gatto [6]. See also [21], in which an anal- 
ogous quantity is used as the invasion fitness of a mutant; a discussion along 
their lines is, however, rather beyond the scope of this paper. 

Comparing the process Z^ with a process Z^ which has bj = bo and 
dj = do , one immediately obtains 

E (1) Z W < e -vt e (b -do-~f)t_ 

Hence, if 

bo — d$ — 1/ <0, 

one has G'(0) < 1. 
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4. Convergence to equilibrium. In this section we prove the convergence 
of the solutions of (1.2) to the unique positive equilibrium, when it exists, 
or otherwise to the extinction equilibrium given by e° = (1, 0,0,... ) T . Con- 
ditions (HI) and (H2) are assumed to hold throughout the section. We be- 
gin with two natural bounds on the mean patch size, the first of which 
bounds s(t) away from infinity. 

Lemma 4.1. Let p° > and let 

oo 

s(t)=J2jPi(t)- 

j=0 

Then 

limsups(t) < +oo. 

t— >oo 

Proof. Multiplying both sides of (1.2) by i and summing for % from 1 
to oo, we obtain 

oo oo 

(4.1) s'(t) = Y^jb jPj {t) -J^jd m (t) - ua(t). 

3=0 j=0 

Note that, as in the previous section, the interchange of derivatives and 
sums is justified, if p° € D(A2), by the fact that the solution p(t) € ^([O, T]; m 1 ) 
and satisfies (3.4). By density, (4.1) then holds for all p° E D{A). 

Now, using the concavity of xb{x) and the convexity of xd(x) as in (3.8), 
we have, from (4.1), 

(4.2) s'(t) < {b(s(t)) - d{s{t)) - u}s(t). 
By standard comparison arguments, we easily obtain 

limsups(t) < s, 

where 

s = inf {s > : b(s) < d(s) + u}. 
The set is not empty because of (H2). □ 

The next lemma gives the complementary comparison result, bound- 
ing s(t) away from when G'(0) > 1 and p° ^ e . Its proof is very much 
more difficult, and is the subject of Section 5. 

Lemma 4.2. Let G'(0) > 1 and d^ < +oo. If p° € C, p° ^ e , then 

liminf s(t) > 0. 

t— »oo 
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Now, if G"(0) > 1, let s* be the unique positive fixed point of G, as in 
Theorem 3.1; if G"(0) < 1, let s* = 0. In the next two lemmas, we show that 
s(t) converges to s*. 

Lemma 4.3. Under the same assumptions as in Lemma 4.1, we have 

limsups(t) < s*. 

t— >oo 

Proof. Assume, if possible, that 

limsups(t) = s > s*. 

t— >oo 

From the proof of Theorem 3.1, we then have G(s) < s. Choose e such that 
G(s + e) < s, and then choose to such that s(t) <s + e for all t > to- 

If we take s(t) as a fixed given function, we see that the solution of (1.2) 
can be interpreted as the distribution of an immigration, birth, death and 
catastrophe process Z(t) with time varying immigration rate s(t), starting 
at time to with distribution p(to)- By an easy stochastic comparison (see [3]), 
that process is dominated in t > to by a process Z^ s+£ ^ with constant immi- 
gration rate s + e and with the same initial condition p(io). In Theorem 2.3, 
it is shown that Z^ s+£ ^ is positive recurrent and that its equilibrium distri- 
bution has finite mean G(s + e) as in (3.2); furthermore, from Theorem 2.3 
and from (2.19) with 5 = 0, it follows that 

(4.3) lim EZ( 5+£ ) (t) = G(s + e) 

t— >oo 

if J2j>ijPj(to) < 00 > true for all p° G T>(A) because of Theorem A. 
Hence, if p(t) is the distribution of Z(t), we have, using also (2.18), 

oo 

lim sup ^2 jpj (t) < G(s + e) < s. 

On the other hand, p(t) is the solution of (1.2) and s(t) was defined as 
J2'jLojPj(t)- The previous inequality thus reads 

oo 

limsups(t) = lim sup ^ jpj (t) < s, 

contradicting s = lim sup^^ s(t). □ 
The companion result is as follows. 

Lemma 4.4. Let G"(0) > 1 and < +oo. If p° eC, p° ^ e , then 

liminf s(t) > s* . 

t— »oo 
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Proof. Assume, if possible, that 

liminf s(t) = s < s*. 

Then, from Lemma 4.2, we have < s < s*. Since G(s) > s, as seen in the 
proof of Theorem 3.1, we can choose e such that G(s — e) > s. 

As in the proof of Lemma 4.3, choosing to such that s(t) > s — e for all 
t>to,we can compare the process with immigration rate s(t) to the process 
with immigration rate s — e. In this way, we obtain 



hminfs(t) = liminf JPj(t) > G(s — e) > s, 

3=1 

reaching a contradiction. □ 

Combining these lemmas, we can prove the following theorem. 

Theorem 4.5. Let (HI) and (H2) hold, and let p° € C\{e }. Then the 
solution of (1.2) converges to the unique positive equilibrium, if G"(0) > 1 
and doo < +oo, and to e° if G'(0) < 1. 

Proof. The previous lemmas together yield 

lim s(t) = s* . 

Now, the interpretation of p(t) as the distribution at time t of an immigra- 
tion, birth, death and catastrophe process Z with immigration rate js(t) 
shows, as in the proof of Lemma 4.3, that p(t) is asymptotically bounded 
between the distributions of the processes Z( s *~ £ ) and Z( s * +£ ) for any e > 0; 
that is, for any / > 0, 

zZ^j ~ £) ^limmf^pjW < limsup^pj(i) ^/Z^j +£) ■ 

3>l ' 3>l 3>l 3>l 

But Theorem 2.4 implies the continuity in s of ^ s \f) with / = 1[{ )00 ), prov- 
ing the theorem. □ 

Remark 4.6. The condition doo < +oo is used in the proof of Lemma 
4.2. There is no reason to suppose that it is necessary for Theorem 4.5 to 
be true, but our proof makes essential use of it. 

5. Repulsion from the extinction equilibrium. The aim of this section 
is to prove Lemma 4.2. To do so, we employ a result from the theory of 
persistence, which we now recall. 

Let X be a metric space (with metric d) which is the disjoint union of 
two sets X\ and X2, and suppose that $ is a continuous semiflow on X\. 
Thieme [25] gives the following definitions: 
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• X 2 is a weak repeller for X\ if 

limsupd(<&t(a;i), X 2 ) > VxiSli. 

t— too 

• X2 is a uniform weak repeller for Xi if there exists e > such that 

limsup(i($ t (xi), X 2 ) > e VxiSli. 

• X 2 is a strong repeller for Xi if 

liminfd($ t (xi),X 2 ) >0 VxiGXi. 

t— >oo 

• A" 2 is a uniform strong repeller for Xi if there exists e > such that 

liminf d(& t (xi), X 2 ) > e VsiGXi. 

In our application, the space X will be the convex set 
C=|pGm 1 :p>0, X^i = 1 j 

with 

00 

d(p,g) = |p - qo\ + ^2j\pj - 

and the continuous semiflow = $>(t,p) is given by the solution p(t) 

of (1.6) with p(0) = p. We take X 2 to be {e } and Xi := C\ {e }; with these 
definitions, the thesis of Lemma 4.2 is that X2 is a strong repeller for X\. 

To prove the lemma, we use Theorem 6.2 of [25], which we state in a form 
simplified to our present needs. 

Theorem B ([25]). Let X be a metric space which is the disjoint union 
of the two sets X\ (open in X) and X 2 ; let & be a continuous semiflow on 
X\. Assume the following: 

(A) There exists a subset Y\ C X\ such that, for all x G X\, there exists 

t(x) > such that G Y\ for all t > t(x). 

{Ce,i) For any y G Y\, the orbit $([0, 00) x {y}) has compact closure. 
(6*6,2) UyeY^iu) h as compact closure, where, as usual, ui(y) is the co-limit 

set. 

(R) The set Y\ PI {x G X; d(x,X 2 ) = e} is bounded. 
Then X 2 is a uniform strong repeller whenever it is a uniform weak repeller. 

We prove that X 2 is a uniform weak repeller, and then Theorem B lets 
us conclude that X 2 is a (uniform) strong repeller, which is the thesis of 
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Lemma 4.2. To start with, we show that the assumptions of Theorem B are 
satisfied. 

Lemma 4.1 shows that, if we choose 

Yi = ^y-Y^m < s|, 

then assumption (A) holds. Indeed, the proof of Lemma 4.1 shows that, if 
y G Yi, then 3>t(y) G Y\ for all t > 0. Assumption (R) is immediate, because 
X2 is bounded. The following lemma establishes the other two assumptions 
of Theorem B. For its proof, note that a set E C m l has compact closure if 
(and only if) limjv— nx>Yl'n=N ^\ x i\ = uniformly for x € E\ that is, if, given 
any e > 0, there exists N = N{e) > 1 such that J2^=n ^\ x i I < e f° r an x £ E. 

Lemma 5.1. If the continuous semi- flow $ is given by the solutions p(t) 
of (1.6) and C, X\ and X2 are the sets defined above, then assumptions 
(6*6,1) and (6*6,2) hold. 

Proof. As in Section 4, observe that p{t) = $i(y) is the distribution of 
an immigration, birth, death and catastrophe process Zt with immigration 
rate s{t) starting at time with distribution y. If y G Yi, this is dominated 
by an immigration, birth, death and catastrophe process Z^ s> > with con- 
stant immigration rate s (because of the previous remark), whose transition 
probabilities we denote by 

p ij (t)=P(Z t ®=j\Z ®=i). 
Stochastic comparison then gives 

OO OO OO OOOO 

(5.1) y n Pn(t) < Y n Yy i P in ^ = Yy*Y n ^n(*)- 

n=N n=N i=0 i=0 n=N 

To estimate the right-hand side, we use (2.15) in Theorem 2.3; choosing 6 
such that c(l + 5) < u, we obtain 

00 -. 00 -. 

Y « <jsY n 1+s PUt) < ^\z^) 1+5 < 

n=N n=N 

uniformly for all t > 0, where Cj is a constant depending only on i. Note 
also that, for 5 = 0, (2.15) implies that 

(5.2) Y n Pon{t) = E®(Z t W) < K t (j + 1). 

n=0 

To prove (C6,i), take y G Yi; choose e > 0. Find N± such that 

i=Ni 1 
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and pick N 2 such that N 2 > (^-) 1 / <5 for i = 0, . . . , iVi; then N 2 is the required 
constant. In fact, using (5.2), we obtain 

oo oo Nl OO OO oo 

J2Vi n Pin(t) =J2Vi np in (t)+ Y Vi Y n Pin(t) 

i=0 n=N 2 i=0 n=N 2 i=Ni n=N 2 

Ni q oo Ni oo 

^X>^r + E y,(i + i)^i<fE^ + 2A 'i E n^ e - 

i=0 2 j=JV a i=0 i=N± 

In order to prove (6*6,2), we prove that, for any e > 0, there exists iV = 
^V(e) > 1 such that, for all y G Yi, there exists to = to(y) such that 



00 00 



Y n E ViPin(t) < e for a11 * > *0- 

n=AT i=0 

Indeed, assume that this is true, and take q G oj(y) for some y G Y\. Then 
there exists a sequence with — ► 00 such that 

00 

(5.3) ^ n|p n (t fe ) - g n | ->• as k -> cx>. 

n=0 

Take such that > io(y) an d that the difference in (5.3) is less than e. 
Then 

00 00 00 

Y n 1n < Y n \Pn{tk) ~q n \ + Y Pn{tk) 
n=N n=N n=N 

00 oo oo 

< Y n \Pri{tk) - q n \ + Y n E^ ire (*) < 2e ' 

n=0 n=N i=0 

using also (5.1), so that (6*6,2) is proved. 

Now choose 5 > such that c(l + 5) < v, and recall as above that, for 
each i, there exists Q < 00 such that limsup i _ f00 E^)(Z t ( s )) 1+<5 < Cj. Hence, 
for each i, there exists i (*) such that E®(Z t &) 1+s < 2Q for all t > t and, 
hence, that 



np in (t) < 2CiN~ & for all t > t {i). 

n=N 

Fix e > 0. Choose ?/ G Yi and find N± = Ni(e, y) such that 

00 

E m < 

i=Ni 1 

where K\ is as in (5.2); set to(y) = maxj=o jvj to(i), an d choose 

N=Ue~ 1 max C< 

l<i<iVi 
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Then, for t > to(y), we have 



X^ E n P™{t) =^2Vi E + E n Pm(0 



i=0 n=N 



i=0 



n=N 



i=N!+l n=N 



i=0 



i=Ni+\ 



proving (5.3). □ 



Now we prove that {e } is a weak repeller through linearization. Since 
we restrict our considerations to vectors p(t) in the convex set C, we have 
Po{t) = 1 — J2j=iPj(t)- Hence, we need only examine the vector (pi,P2, ■ ■ ■ ) T ■ 
With a slight abuse of notation, we now set 



< +oo 



X= x=(x 1 ,x 2 ,...) T €£ 1 :^j| 
I j=i 

with norm ||x|| = J2j=ij\ x j\i noting that e° now translates into the point 
of X, and we denote here by A and F the operators defined in (1.4)-(1.5) but 
restricted to X, and using po = 1 — J2j=iPj i n t ne definition of F. We then 
define X + to be the nonnegative cone in X; note that X + is the counterpart 
of the convex set C defined above. 
Equation 



(5.4) 



p = Ap + F(p) 



corresponding to (1.6) now has as the equilibrium, corresponding to the 
extinction equilibrium e° of (1.2). We again use ^t{u°) to denote the solution 
of (5.4) satisfying u(0) = u°. This corresponds to the semi-flow of Lemma 5.1, 
except that we now neglect the Oth component. Note that the metric in the 
convex set C is equivalent to the norm in X , since 



d(u, v) 



u 



oo 

■vq\ + ^2i\ui 
1=1 



i=l 
oo 

< 2^2i\m - 

i=l 



i=l 



oo 

i=i 



Vi 



2\\u - v\\ x , 



while obviously ||tt — v\\x < d(u,v). 

Note also that A is the generator of a defective Markov process, the pro- 
cess of Section 2 restricted to the state space N \ {0}. In the rest of 
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this section, we only consider processes with zero immigration rate; thus, 
when there is no ambiguity, we drop the superscript and denote by Z% 
the process with zero immigration rate. 

From the results of Section 2, one immediately sees that Zt is exponen- 
tially absorbed at 0; more precisely, (2.13) with f = e and I = s = implies, 
in the present notation, that 

(5.5) \\e At \\ <Ce~ at 

for some positive constants C and a. This implies that 

{IRA > -a} C p(A); 
moreover, we have the representation 

(5.6) ((A-^)~H = E^( A )' 

3 

where "~" denotes the Laplace transform and Pji(t) is ¥(Zt = i\Zq = j). 

We now discuss the stability of the equilibrium of (5.4) using the lin- 
earization principle. We first note that 

(5.7) F'(0)u = ip(u)e 1 VtiGl, 
where 

(5.8) (p{u)=~ 1 ^ju j 

j 

and e 1 = (1, 0, 0, . . . ) T . Since ^'(0) is one-dimensional, hence, compact, the 
essential spectrum [27] of A + F'(0) coincides with that of A, which, from 
(5.5), is less or equal than —a. The type of the semigroup e^ A+F can 
then be understood from the spectrum of A + -F'(O). 

Using (5.7), we can establish, through direct computation, the following 
result. 

Lemma 5.2. If \ is in p(A), then A belongs to p(A + F'(0)) if and only 
if ip{{\ — A)~ 1 e 1 ) ^ 1. In that case, 

(5.9) (A - A - F'(0)y\ = (A - A)~ l v + ^((X-A)-^ ~ A) ~ V - 
On the other hand, if 

(5.10) ^((A-A)-V) = l, 

then A is an eigenvalue with corresponding eigenvector v = (A — A)~ 1 e 1 . 
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From this lemma, we see that an important role is played by the roots 
of (5.10) in the half-plane {MA > —a}. Using the representation (5.6) and 
standard results on the Laplace transform, as used, for instance, in the 
theory of age-dependent populations [16], we have the following: 

Lemma 5.3. There exists, at most, one real root Xq > —a of (5.10). If 
Xq exists, all the other roots X satisfy MX < Ao,' if there is no real root, there 
are no complex roots in {MA > —a}. In any strip {a < MA < b}, there are, 
at most, finitely many roots. 

Finally, if 



then Ao > [=] 0; on the other hand, if Ro < 1, if there is a real root Ao, it 
satisfies Ao < 0. 

Remark 5.4. Note that 



with G as given in (3.2). 

From here on we assume that Ro > 1. Hence, the real eigenvalue Ao is 
positive. We denote by Ai, . . . , A& (with k > 0) the other roots of (5.10) such 
that MXj > 0, and by Afe+i, ■ ■ ■ ,X n (with n > k) the roots such that MXj = 0. 
Since the continuous spectrum (if it exists) of A + F'(0) is contained in 
{MA < —a}, we can split the spectrum of A + F'(0) in three spectral sets 



a u = {Ao, Ai, . . . , A fc }, <? c = {Afc+i, . . . , A n } and a s = {A € a(A + F'(0)) :MA < 



By standard results (see Theorem III. 6. 17 in [17]), X can be split into the 
direct sum of three subspaces X u , X c and X s , all invariant under A + F'(0). 
Moreover, X u and X c are finite-dimensional {X u includes at least vq, the 
eigenvector corresponding to Ao, while X c may well consist only of 0). This 
would be enough to establish the instability of the equilibrium. However, 
we wish to prove that all initial data u > 0, u^O, are repelled away from 0, 
and this requires further work. 

The following lemma uses the results of Bates and Jones [4] to establish 
the existence of unstable and centre stable manifolds W u and W cs for equa- 
tion (5.4) at 0. The conditions of their Theorem 1.2 are satisfied in view of 
Arrigoni's results, as summarized at the end of Section 1, together with (5.5) 
and the properties of the eigenspaces discussed following Lemma 5.3. 

Defining X cs = X C (B X s , and letting P u and P cs denote the corresponding 
projections, [4], Theorem 1.2 and its consequence (P3) yield the following 
result. 





0}. 
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Lemma 5.5. There exist a neighborhood U 3 and Lipschitz functions 
h u . p U( jj^ _^ X cs and h cs . pcsQjj X u with h u( Q j = (/j«)'(o) = h cs (0) = 

(h cs )'(0) = such that 

W u = {u u + h u (u u ) : u u 6 P U {U)} 
is the unstable manifold (in U) ofO, and 

W cs = | u c S + h cs^ u csy u cs £ pcs^j 

is a centre-stable manifold. 

Furthermore, there exists a neighborhood V C U of such that, if u° G 
V \ W cs , then there exists r > such that <£ T (u°) ^ V. 

The final statement of the lemma shows that, if a solution comes close 
enough to to be in the neighborhood V, and if it is then at a point not 
in W cs , then it has to leave V at some later time. Hence, the limes superior 
of any solution curve is necessarily positive, if it can be established that, for 
some e > 0, no points of X + n B e except for are in W cs , where B £ denotes 
the ball of radius e centred at 0. If this is the case, then {0} is a uniform 
weak repeller for X + \ {0} in the system (5.4), which is equivalent to {e } 
being a weak repeller for C \ {e } in (1.6). Applying Theorem B, Lemma 4.2 
would then follow. 

To show that indeed W cs D X + n B £ = {0} for some e > 0, we begin by 
writing the eigenprojections explicitly. 

Lemma 5.6. The projection Pq on the eigenspace corresponding to Xq is 

y((\ -A)- l v) v 
Pqv TTT\ ^rT^ A ° -A ' ) e - 

The projection P u on X u is given by 

where Tj is a circle around Xj that does not include other elements of the 
spectrum. 

Proof. It follows from the construction of the projection operators as 
in formula (III. 6. 19) of [17] and from (5.9). □ 



Note that 

^((Ao-A)-V; 



POO °° 

/ te-^YiPui^dtKO. 
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On the other hand, it may well be <p'((\j - A)~ l e l ) = when l<j< n, so 
that the other projections may have a more complex form. 
As a consequence, we immediately have the following result. 

Lemma 5.7. IfveX cs , then <p{(X - A^v ) = 0. 

Proof. The explicit representation of Po shows that if ip((\o — A)~ l v) ^ 
0, then Pqv ^ 0. However, v G X cs implies that Pqv = 0. □ 

This lemma implies that X cs DX + = {0}, because if v > and v ^ 0, then 
(A - A)~ 1 v > by (5.6); hence, it follows from (5.8) that <^((A - A)~ 1 v) > 
also. This is almost what we need, since W cs and X cs are close to one another 
near 0, and we are thus close to showing that W cs n X + DB £ = {0} for some 
e > 0. To make the transition from X cs to W cs , we first show that, for v > 0, 
<^((Ao — A)~ l v ) is large enough. 

Lemma 5.8. Assume that d^ < oo, and take v > 0. Then 

(5.11) ^((A - A)-M > - M— - . 

«oo + 7 + 17 + ^0 

Proof. We start from the identity 

/■OO 

V^(A )= / e- Ao *E^(Z t )cft. 

An easy coupling argument shows that is stochastically larger than a 
death-and-catastrophe process with death rate d^ + 7. Hence, 

rOO 



/ e- x °%(Z t )dt>j 
Jo Jo 

Now, if v > 0, we have 

<^((A - =^2iv j P ji (Xo 



(A +7+rfoo+^)t ^ _ J 

A + 7 + doo + ^ ' 



> 



.7 



A + 7 + (i oo + i/ doo + 7 + z^ + Ao □ 



Using the above lemma together with Lemma 5.7, we can now show that 
the norm of v~ is quite large, whenever v € X cs . Here, v~ denotes the nega- 
tive part of v : v = v + — v ~ , with (v + )i = max{0, Vi} and (v~)i = max{0, —Vi}. 

Lemma 5.9. Assume that d^ < 00. // v € X cs , there exists n > such 
that \\v~ || > i]\\v\\ . 
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Proof. From Lemma 5.7, we have 

= ¥>((A - A)~ l v) = V9((A - A)'^) - p((Ao - A)~ x v~) 

using Lemma 5.8 and the obvious identity \\ip\\ = 1. Hence, using ||tf 
\\v\\ — \\v~\\, 

(A " A)" 1 !) + - ; 1 —) \\V~\\ > 



which yields the thesis. □ 

We now use this result, together with the closeness of X cs and W cs , to 
conclude that W cs nX + nB £ = {0} for some e > 0. 

Lemma 5.10. Assume that doo < oo. Then there exists e > such that 
v>0, v G B e (~l W cs implies v = 0. 

Proof. First take 5 such that ||t> cs || < 5 implies ||fo cs (i> cs )|| < §||-u cs ||. 
Then take e = 5/||P cs ||. Assume that v = v cs + h cs (v cs ) > with ||«|| < e. 
Then it follows that ||n cs || = \\P cs {v)\\ < S. 

Split v cs = (v cs ) + — (v cs )~ . Then we have 



i:u? s <0 i:«f<0 i=l i=l 

= -IIKTII + ll^ c > cs )ll < -^Ibl + §K1, 

using Lemma 5.8 and ||u cs || < 5. This contradicts with u > unless v cs = v = 
0. 

□ 

We have now proved what we need to show that {0} is a uniform weak 
repeller for X + \ {0}. The details are as follows. We recall that we have 
R = G'(0) > I. 

Lemma 5.11. Assume that d^ < oo. Then there exists £q such that for 
all u° >0, uV 0, 

limsup || $t (it ) || > so- 

t—*oo 
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Proof. Take 

£n = Amm<£, inf \\v\\>, 
Z { v£X\V J 

where e is that of Lemma 5.10, while V is that of Lemma 5.5. 

Assume that )|| < 2sq for all t > to. Since u° > 0, the invariance of 

the positive cone under (5.4) gives moreover, $ t (u°) Hence, 

Lemma 5.10 implies that $4 (u°) ^ VF CS . From Lemma 5.5, it then follows 
that $ T (u°) £ V for some r > t , contradicting ||$ t (u°)|| < 2e for all t > t . 
□ 

PROOF of Lemma 4.2. Going back to the semi-flow $ t on C, note that 

oo oo oo 

d($ t (u ),e°) = |1 -po(*)| +J2i\Pi(t)\ =^Pi(t) + £«Pi(*) 

j=i j=i j=i 

oo 

<2^ip i (t) = ||* t (uo)|| I 
i=i 

while obviously 

||*t(«o)||<d(*t(«o),e°). 

Hence, Lemma 5.11 states that {e } is a uniform weak repeller for C\{e }. 
But now Theorem B, together with Lemma 5.1, yields Lemma 4.2. □ 
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